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We study master variables in the Regge-Wheeler-Zerilli formalism. We show that a specific choice 
of new variables is suitable for studying perturbation theory from the viewpoint of radiation reaction 
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I. INTRODUCTION 

Black hole perturbation is a powerful tool for the 
evaluation of gravitational waves from a binary sys¬ 
tem when its mass ratio is large 0 1 i Although 
any systematic method to calculate the radiation re¬ 
action to the particle motion has not been established 
so far, there are various new developments in this field 



mal prescription to extract the self-force was developed 
in [4j, generalizing the work of DeWitt and Brehrne [jj 
on electromagnetic self-force to include the gravitational 
case. These results were further verified by an indepen¬ 
dent, and different, axiomatic approach by Quinn and 
Wald [{J. The prescription can be summarized as fol¬ 
lows. The retarded held can be obtained in terms of 
Greens functions which can be formally decomposed into 
“direct and “tail parts. Roughly speaking, the “direct 
part of the held is that part which has support only on 
the future light cone, emanating from the source point. 
The “tail” part is composed of contribution due to cur¬ 
vature scattering which pervades inside the future light 
cone of the source point. The analysis presented in 0[3 
indicates that, the particle motion, after taking into ac¬ 
count the self-force, follows a geodesic on the geometry 
perturbed by adding the “tail” part to the original back¬ 
ground spacetime. 

The actual isolation of the “tail” part is not an easy 
task. There are ways to calculate the “full” Greens func¬ 
tion but there is no direct method to compute the “tail” 
part alone. Hence, the standard prescription that has 
emerged in the past few years is to subtract the “di¬ 
rect” part from the “full” metric perturbation. Here 
lies the well known “gauge problem”. In the standard 
methods for constructing full metric perturbation, we 
first solve the equation for master variables and then 
from these master variables we reconstruct the metric 


perturbations. The result is naturally written in a spe¬ 
cific gauge such as Regge-Wheeler (RW) 19|, 20] or radi¬ 
ation gauge HmHEl. ^ ie °th er hand, the “direct” 
part is evaluated in the harmonic gauge associated with 
the particle trajectory. Therefore, before any meaning¬ 
ful subtraction we need to relate these expressions which 
are in different gauges. This is by no means an easy task 
since we do not know the necessary gauge transforma¬ 
tion a priori. This additional task to find the appropri¬ 
ate choice of the gauge parameters makes the problem 
much harder to solve; this is the aforementioned “gauge 
problem”. The attempts for subtraction of the “direct” 
part in the RW gauge were reported by Mino Mino 3 
and Sago et. al. jig . 

In this paper we would like to revisit the problem 
of metric perturbations reconstruction from the master 
variables in case of Schwarzschild background. In this 
approach, based on Regge-Wheeler-Zerilli formalism, the 
problem of metric reconstruction is relatively well under¬ 
stood 3E3. What we would like to discuss here are 
possibilities of improvements on this formalism. 

It is well known that in the Schwarzschild case the odd 
and the even parity perturbations naturally decouple due 
to spherical symmetry. Assuming a time dependence of 
the form exp(— iuit), the perturbations of a Schwarzschild 
black hole can be described by a master equation, for 
each partial wave mode decomposed in terms of spherical 
harmonics, as 

^ + (“ 2 -v)<; = s(t, v ). ( i ) 

Here £ is the master variable, S is the source term com¬ 
posed of the matter energy momentum tensor T^ v and 
r* = r + 2M ln(r/2M — 1) is the usual tortoise radial co¬ 
ordinate. The metric components in the RW gauge, h RW , 
are obtained by applying certain differential operators on 
the master variable and on the energy-momentum tensor 
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as 

h RW =U m \C) + U t \T^). ( 2 ) 

Then, the formulas for the metric reconstruction in the 
scheme presented in the original papers [11113 contain to 
in the denominator of the expressions for /d T ). Although 
u) is just a number in the frequency domain, if in denom¬ 
inator it can be an obstacle in computing the metric in 
the vicinity of a particle orbiting a black hole. Suppose 
that the particle moves between r m j n and r max . The ap¬ 
pearance of io in the denominator means that h^ T ' ) is no 
longer localized on the radial shell where the particle or¬ 
bit lies. Instead, the source is distributed continuously in 
the region between r m in and r max . Therefore, the metric 
components are not completely determined by the notion 
of the master variables in this region even if the concerned 
field points are off the shell. In the computation of the 
self-force, the gravitational held exactly on this shell is 
unnecessary. A limiting value evaluated along, e.g., the 
outer radial direction is sufficient for the purpose of com¬ 
puting the self-force. If we can modify the formulation 
so that h( T) is localized on the shell, then we can ap¬ 
ply the formula for the metric reconstruction outside the 
source, which is much simpler. For the even parity case, 
an improved master variable has already been introduced 
by Moncrief ( .24], see also |2flj [ 3 ( 1 ] ). We give here the 
general metric reconstruction formulas in the presence of 
sources, which have not been given explicitly yet, as far 
as we know. We show that u can be removed from the 
denominator by using Moncriefs master variable. Same 
argument follows for the odd parity case, i.e., by intro¬ 
ducing a new improved master variable, we can remove 
appearance of u> in the denominator. Complete expres¬ 
sions for the metric reconstruction are also presented for 
this case. 

Another complication which arises is from the well 
known fact that the potentials for odd and even parity 
cases differ from each other. The potential for the even 
parity case is, relatively, much more complicated. Hence, 
it would be useful if we could formulate the even parity 
perturbations to satisfy the same master equation with 
the odd parity case. Chandrasekhar has already given 
a unified approach, known as Chandrasekhar transfor¬ 
mations, and shown the relation between RW and Zerilli 
equations [26| (for a comprehensive review see [2?}). In 
this paper we also derive the full metric reconstruction 
formulas for the even parity perturbation by using the 
master variable obtained via the Chandrasekhar trans¬ 
formation. Under the requirement for this new master 
variable to satisfy the RW equation in vacuum, we can 
still modify its definition by adding a combination of the 
metric components which appear on the left hand side of 
Einstein equations, since it is zero in vacuum. Examining 
all the possibilities of such a modification, we have con¬ 
cluded that we cannot eliminate u> from the denominator 
in the expression for h^ T \ Unfortunately, as it turns 
out, no dramatic simplification happens by reformulat¬ 
ing the formulas solely in terms of the variable obtained 


via Chandrasekhar transformation, although the impor¬ 
tance of this transformation is not reduced at all by this 
fact. 

The paper is organized as follows. In Sec. II we discuss 
the equations for both odd and even parity cases. We 
have provided the explicit expressions for source terms 
corresponding to the new master variables. In Sec. Ill 
the even parity master variable which satisfies odd par¬ 
ity homogeneous master equation is discussed. We briefly 
summarize the results obtained in this paper in Sec. IV, 
with a speculation towards an alternative method to com¬ 
pute the regularized self-force subtracting the direct part 
at the level of the master variables. 


II. IMPROVED MASTER VARIABLES 


We begin with reexamining the Regge-Wheeler-Zerilli 
(RW) formulation. In this formalism a master equation 
for a master variable is derived, which are called the RW 
equation and the RW variable, respectively. Once we 
know the solution for the RW variable, all the compo¬ 
nents of the metric perturbation can be derived from it. 
This scheme is well known. What we have shown here 
is that it can be improved, in the sense discussed earlier, 
by introducing alternative master variables. 

We consider the Schwarzschild metric, 


ds 2 




+ r 2 (d,0 2 + sin 2 9 dip 2 ) , 


( 3 ) 


as the background. The 10 metric components can be de¬ 
composed into “odd” and “even” parity modes. We use 
the notation in which, after harmonic decomposition for 
the angular dependence, H 0 , Hi, H 2 , h 0 ,h\ , K and G 
are the components of metric perturbations for the even 
parity modes, and ho, hi and h 2 are for the odd par¬ 
ity modes. Here we assume that the time dependence is 
given by exp(— iuit). Similarly, the components of the en¬ 
ergy momentum tensor can be decomposed into odd and 
even parity modes. A^\ A, B®, B, and F are 
the expansion coefficients for the even parity modes, and 
Q^°\Q and D are for the odd parity modes (we follow 
throughout notation of Zerilli for the metric perturbation 
and the e nerg y-momentum tensor with slight modifica¬ 
tions; see [Uj for the basic equations such as the law of 
gauge transformation and the definitions of the harmonic 
expansion coefficients of the energy momentum tensor). 


A. Odd parity 

First, we consider the odd parity case. The RW gauge 
choice corresponds to setting h^ w = 0. Here, the vari¬ 
ables with a superscript RW means the quantities are in 
the RW gauge. The nontrivial set of Einstein equations 
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for the odd parity mode is 
t rw , * w /.2 RW 


"0,rr + ^( r h l ),r + 


AM 

-2(1 +A 

r 


< w 


r(r — 2M) 


8tt 


Vl + X r - 2 M 


Q (0) . 


( 4 ) 


h RW 9 

+ iu}h r + 2A ( r - - rh™ 


87 ri 

v/TTa 


(r - 2M)Q , 


( 5 ) 


1 - 


2M\ 


h 


r J 


RW 

l,r 


+ iw [ 1 — 




-1 


l fltv , 244 . itw 
n o + „2 “1 


87TJ 


■\/2A(l + A) 


r 2 D. 


( 6 ) 


^From the above equations, and using the conventional 
gauge invariant master variable ■“•’y 


For reconstruction is straight forward since there is 
only m M \ In reconstructed h RW first term in the square 
brackets corresponds to and second term is h ( - T \ 
The trouble with the expression for h^ is the presence 
of u> in denominator, as anticipated earlier. Even if we 
rewrite this expression using conservation law 

V2A D = + (3 - —) Q + (r— 2M)Q r , 

r — 2 M \ r J 

( 13 ) 

this uj cannot be removed. This fact implies that we need 
time integration of the source term in the reconstruction 
of metric perturbation. Hence the reconstructed metric is 
not solely determined by the master variable even when 
the energy-momentum tensor vanishes on the spherical 
shell containing a given field point. 

We, therefore, introduce a new gauge invariant variable 
defined by 


(o). 


r — 2M 


X = 


K 


RW 


( 7 ) 



—iujh 


RW 

1 


- h RW 

n 0,r 


+ h 

r 


RW 

0 


(14) 


we can derive a second order differential equation as 

[t d 2 . +oj 2 - V RW {r)\ (°\ = S^ x . (8) 

This is the well known Regge-Wheeler equation |r9|. Here 

V™ = (1 - “) (r 1 - “) ' < 9 > 

is the Regge-Wheeler potential and the source term is 
given by 


Using the definition Q with one of the odd-parity field 
equations, we can verify that 


- »+ = + + 8 " r(r +H W . (15) 

2AVTTA 


Hence, is equivalent to time integral of the original 
variable outside the source distribution. The expres¬ 
sion applicable to an arbitrary gauge has the same func¬ 
tional form as in the RW gauge: 


£>+ 87TI 

5 “ VXTT 






Q 


( 10 ) 


Here A is defined in terms of t, the total angular momen¬ 
tum of the spherical harmonics, as 


(£-!)(£ + 2) 


( 11 ) 


Once the master variable and the energy-momentum ten¬ 
sor are given, we can reconstruct the metric perturba¬ 
tions. To distinguish these reconstructed variables from 
the original ones, we associate them with an over-hat in 
the same way as and hS T \ Combining the Einstein 
equations using the definition of the master variable, the 
necessary formulas for the reconstruction can be derived 
as 


h RW 

h RW 

rl 0 


»- 


r-2 M 


1 - 


X, 

“)[ 

8nir 2 ^ 

4— D 

\/2A(A + 1) 


1 

iu) 


(' r (o) x), r 


( 12 ) 


io) c= _l_ 
^ 2A 


—iuihi — ho, r H —ho 


Substituting Eq. (□ into Eq. © , we recover 

[d 2 r ,+u 2 -V RW (r)\ r = S (0) ^ 
with new source term 

_ 87r(r - 2 M) 


2AV1 + A 


wrQ — d r (rQ^) 


(16) 


(17) 


(18) 


Here we have used conservation law m to simplify the 
expression. The source term S ( ^ does not have time 
integral although <0 \ is a time integral of the original 
variable This is expected a priori. If the source term 
for has integration constant then it is not uniquely 
determined, which contradicts with the fact that it is 
a gauge invariant variable. Here, in order to illustrate 
the way how we found the new variable, we took rather 
lengthy steps to obtain the master equation ra passing 
through the equation for the original master variable ©• 
But, of course, one can directly verify the final result by 
combining first two odd-parity equations © and ©• 
Now we consider the reconstruction of the metric com¬ 
ponents from this master variable There are two 
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non-vanishing components h RW and h RVi in the RW 
gauge. They are to be solely determined from (°'£, if the 
metric perturbation satisfies Einstein equations. From 
Eq. J3J) and the definition of we immediately have 


h RW 


icor 2 
r - 2 M 


(o) C + 


Airir 3 ^ 

\VT+x Q ' 


(19) 


Once we know h RW , we can reconstruct h RW by using 
Eq. 0 as 


h^ w = (r—2M) + i (o) C + a ^ Q (0) ) • ( 2 °) 


This time, the w T ) part does not have to in the denomi¬ 
nator. Therefore, one can simply set the source terms to 
zero to obtain the formulas for the reconstruction of the 
metric perturbation in vacuum region. We notice here 
that h( M ' ) is also free from annoying factor w. These 
two facts are actually related. By definition, the defin¬ 
ing expression for a gauge invariant master variable does 
not have u> in the denominator. Otherwise, the gauge 
invariant variable would be ambiguous due to integra¬ 
tion constant, and information of metric perturbations 
in the vicinity of a spherical shell, specified by t and r, 


will be insufficient to determine the corresponding gauge 
invariant variable there. Let us assume that also 

does not have u> in the denominator. In the vacuum 
case, we can consider a cycle of operations starting with 
h, going through the master variable, and again coming 
back to h by using h^ M \ Throughout this cycle, there 
is no to in the denominator. Hence, if a homogeneous 
solution of metric perturbations including its derivatives 
near a spherical shell is given, this cycle should reproduce 
the original metric perturbations. The formulas compos¬ 
ing this cycle will not change even if there exists matter 
source away from the spherical shell. If h < ' T 1 has to in 
the denominator, this term gives an additional contribu¬ 
tion even if matter source does not exist there. This is 
a contradiction. Hence, if does not have u> in the 
denominator, neither does hi T \ 


B. even parity 

Next, we look at the even parity case. The RW 
gauge choice corresponds to setting h^ RW = h[ e ^ RW = 
qRW _ q The set of field equations for the even parity 
modes, with the source terms, is 


J 


1 - 


2 M N 

r J . 


1 - 


2 M 


(K R ™ - 1 Kr)+ (3 — 


RW\ 


5 M\ 1 


-K 


RW 


- 4 ( H RW -K 


■~RW \ 


= -8ttH (0) , (21) 


-iwK RW - iuj—{K RW - H rw ) + i* , 

' r r v 2 ’ r(r-2M) 


~ (H RW + K RW ) 

r z 

k rw __ 0_±V H RW = _^ A d) j (22) 

r 2 v 2 


(r - 2 M) 


—u> 


(r - 2 M) 


K 


RW 


- (l - K RW + 2 iwH RW + (r 

, (1 + A) 


RW _ X RW ) 




RW 


(.k hw - h rw ) 


+ iu(H HW + K hw ) = 


= -8? tA, (23) 
8ni 


Vl + X 


rH (0) , (24) 


2 M 


u 


r 

-l 


2 M 


iu>H RW +1- (H RW - I< KW ),r + — 


tRW 


1 


M 


- i - — (Hr - k w ) = 


8n 


VT+X 


(r - 2 M)B, (25) 


1 - — ) (K rw + H rw ) + 


r) 


t rw 


2(r — M) RW _ 1 
r(r-2M) 1 r 


1- 


M \ 


r J 


1 H?r - 


[I<r - H™] + (i - y) \k rw - 2 icoH R 

^1 + Hr - ^+ 2 X \ h rw - H rw ) = 8\/2t rG (s) , (26) 

r 2 F. (27) 


ttRW ttRW _ 

H 0 ^2 — 


167T 


\/2A(l + A) 


The original Zerilli’s master variable, Rim, is defined by 


Rlm = - (e) X=~ 

U! UO 


r — 2M 
X r + 3M 


icor 


r — 2M 


-K 


RW 




RW 


( 28 ) 


and has an ambiguity due to an integration constant. 
We will work, instead, with its gauge invariant form 
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x ■ With the field equations above it obeys the wave 
equation *2£| 


[d 2 r , +uj 2 - V z (r)]( e) X = 5 (<e)x) . 


(29) 


Here, 


f 2 M\ 

v w = y 1 - —) 

2 A 2 (A + 1 )r 3 + 6A 2 Mr 2 + 18A M 2 r + 18 M 3 


r 3 (rA + 3M) 2 

is the Zerilli potential and the source term takes the form, 
(r - 2M) 2 


S( (e) x) = 


_ = H (0) 

(rA + 3M)v / lTA ,r 
(r - 2M)(—12M 2 + 9Mr + r 2 A) (0) 

rV 1 + A(rA + 3M ) 2 

(rA + 3M) 2 
(r - 2M) 2 


-A + 


(rA + 3M) ‘ 
(r + 2M) 2 


(rA + 3 M) (rA + 3M)vTTA 
i/A(l + A) 


B 


(31) 


The formulas for the metric reconstruction are derived 
by combining the Einstein equations using the definition 
of the master variable. Since this is a known result, we 
just quote here the explicit reconstruction formula for 
K rw as an example: 


fcRW _ 1 
L0 


- l- 


2 M 


(e). 


X,r 


r 2 A + (rA + 3M)(rA + 2 M) 


r(r — 2 M) 


r 2 (rA + 3 M) 


(e). 


X 


w(rA + 3M) \ 


= A« 


1 


vTTa 


H (0) 


(32) 


As in the odd parity case, the first term is and the 
second term is F T \ The presence of to in denominator 
in the expression for h <Tl is a signal that this ^x is not 
the most convenient choice of the master variable. 

Analogous to the odd parity case we now define a new 
time integrated variable using vacuum field equations as 


(e)^ = 


r(r — 2 M) 


(A + l)(Ar + 3M) 


[H^ w - rK J 


RW 


+ ■ 


rA + 3M 
r - 2 M 


K 


RW 


.(33) 


In fact, the same variable has been introduced earlier by 
Moncrief 0 (See also Gleiser et. al. j2f|). It can be 
easily checked that satisfies a similar wave equation 


[d 2 , +w 2 -y z (r)] (e) C = S c , 

with a modified source term 


(34) 


^ A (0) - w-^=A (1) + (i + A)(r - 2 M)A 
+VTTA(r - 2 M)B - \j 2i ' l + A) (rA + 3M)F .(35) 


rA + 9 M 
rA + 3 M 


Here, for simplification, we have used the three constraint 
equations, corresponding to T^ ;At = 0, which are 


(30) A (i) = 


1 


V2wr 2 


A r — 


B r = 


r — 2M 
+\/2(l + A )B (0) 

1 


i--!#) 
r — 2 M V r 


M 


(36) 


- 


Mr 


r-2M L\/2(r- 2M)“ (r - 2M)' 

+ M-2r A + + V2G (s) 


r A(°) 


(37) 


1 


r-2M 

-V2(1 + A)G C 


4M 


-^—B^ _ 3 - — B + V2AE 

r — 2Af V r 


(38) 


As explained in the odd parity case, the source term for a 
gauge invariant variable does not have u> in denominator. 
Now we come to the reconstruction of the metric com¬ 
ponents using this new master variable in the RW 
gauge. There are four nonvanishing components in the 
even parity case, namely, K RW , H RW , H RW and H RW . 
We can rewrite them in terms of the gauge invariant vari¬ 
able (e) C as 


RW 
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jjRW 


A(A + l)r 2 + 3AMr + 6 M 2 (e)/ _ | / _ 2 M \ ^ 
r 2 (rA + 3M) ^ V r j U 
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( 39 ) 
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(1 + A)(rA + 3M) 



/ rA + 2M 
V r - 2M 


These reconstruction formulas are local and do not re¬ 
quire any time integrations. 
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III. EVEN PARITY MASTER VARIABLE VIA 
CHANDRASEKHAR TRANSFORMATION 


In this section we have examined the even parity mas¬ 
ter variable that satisfies the same homogeneous mas¬ 
ter equation as the odd parity one. The method to ob¬ 
tain such an even parity master variable is well know 
as Chandrasekhar transformation [2(|. Here, we give a 
short derivation of this transformation, and discuss the 
metric reconstruction scheme using this new variable ^ e \. 
As mentioned earlier, a part of motivation is the useful¬ 
ness of master variables which satisfies the same master 
equation for both the parities. In particular, the master 
equation is much simpler for odd parity case. Another 
point is the appearance of the factor l/(Ar + 3 M) in the 
RW potential in the even parity case, which is absent in 
the odd parity case. This factor mathematically means 
the existence of a singularity at r = —3M/A in the master 
equation. However, this singularity will not be a physical 
one because of the symmetry between even and odd par¬ 
ity cases. This factor l/(Ar + 3 M) is inherited in many 
places of the whole reconstruction scheme. Although not 
a serious obstacle in actual computation, we can expect 
that the reconstruction scheme might simplify a lot by 
using the new variable ( e %. 

Our quick derivation of is based on the 

fact that Weyl scalar contracted with null tetrad 
pb = -Cabcdl a m b l c m d satisfies the same homogeneous 
equation irrespective of the parity p^|. Here l a 
and m° are outgoing and angular null tetrad vec¬ 
tors, respectively. For explicit calculations, we use 
( l a ) = (r — 2M) _1 (r,r-2M,0,0), and ( m a ) = 

(Vorsin#) -1 (0,0,sind,z). The following formulas are 
obtained just by plugging in the explicit metric form into 
the definition of ip, 




and 


2r 2 (r — 2M) 


[(r - 2 M)h$J + rh™ - icorh f 


•RW 


(2 M + iur 2 ) uRW 
r — 2M n ° 


(40) 




(H* w + H* w ) 
r(r — 2 M) 


(41) 


J 


Here the angular dependence, which is given by the spin 
weighted spherical harmonics, is suppressed for brevity. 
We use the same notation ^'Up to represent the coeffi¬ 
cients of Fourier harmonic decomposition, but it will not 
cause any confusion. 

Substituting Eqs. con and m , we can rewrite 
Eq. m i n vacuum as 

(0 V = 

+(3 M - (A + 1 )r)(r - 2M)}^°\(r) 

+r(icur 2 + 3 M — r)(r — 2M)^*( r (r)]. (42) 
Here we have used the field equations for simplification. 
From the equation above and with the aid of Eq. CGJ in 
vacuum, we can express the master variable in terms 
of and its derivative as, 


(o) C( r ) = C[ ( °V] = 


[{urr d — 5 iuiMr 


2(3iu>M + A(A + 1)) 

—AM + r + A(2 M — — (r — 2 M)(iujr 2 

+3 M — r)(°\p ir (r)\. (43) 


Then, with an arbitrary constant C, 

(e) C = cc[ (e V] 


(44) 


should satisfy the RW equation, i.e., the same equation 
that satisfies except for the source term. After a 
straightforward calculation, we obtain 

(e) c>) = 2(r—2M) ( H? w - rK« w + -^L—K RW ^j , 

(45) 

with 

C = 4(3iujM + A(1 + A)), (46) 

which is known as the Starobinsky constant[27^. It is 
also easy to check directly that this new master variable 
satisfies the homogeneous RW equation with the same V 
as the usual RW potential. 

In general case with the source term, we have 


[d 2 r ,+u 2 -V RW {r)]^ = S^, 


(47) 


with 


= (r - 2 M) 


-2 rAW + 2 M - r(1 + A) A (0) + y/7urrAV + 2 [6M - r(l + A)](r ~ 2Af) ^ 
’ i — 2 M r- 


+2 


\6M — r(l + A)](r — 2M) 
y/l + A r 2 


B - 2\/2 


[6 M 2 — Ar 2 (l + A)] 
\/(l + A)A r 2 


\/(l + A)Ar 


• (48) 
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The metric reconstruction formulas for the are given by 
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—r( 1 + A)(rA + 3 M)A (0) + -VlTXM(r - 2 M) 2 {B + VlTXA} - -^=wrM(r - 2 M)A (1) 
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16 [ Vr 2 (rX + 3M) , n1 [3Mw 2 r + A 2 (A + l)]r 2 m 5 

icFh" (r-2M) A()+sl -yr—-• 4< ’ +3, “ MP(r 

_ .M^(rA + 3M), .J,^ + 3 MP( r -2M)|,.,- (rA + 3M) P(e ,- 

^(1 +A)A 2r 2 (r — 2M) s 2r 


} 


ff 2 = 


16 

jcj 2 [(r — 2M) 


2 C> A r 2 V 


A^+i 


V2 


A (1) - A(r - 2Af)P{(l + X)A + VTTXB} + \/2(l + X)XV(rX + 3 M)F 


[O — X(m 2 r 3 + M(1 + X))\P ( e )~ _ O (e) ~ 

2(r — 2M)r 2 ^ 2r 


//o = // 2 


16-7T 


\/2A(l + A) 


r z F.. 


(49) 


Here, T 5 = 3M — r(l + A) and O = 3 Muj 2 r 2 + A(A + 
1)(3 M — r). If we are working in frequency domain 
only, the above choice of master variable is not a bad 
one because of the common potential in master equa¬ 
tion. Whereas, in the time domain we will need time 
integrations for the metric reconstruction due to the fac¬ 
tor |C|- 2 . 

We can modify the master variable by adding combi¬ 
nations of metric components which appears on the left 
hand side of the Einstein equations. Let’s denote these 
combinations by G (.-a so that the Einstein equations are 
formally written as Gu\ = TVq (* = 1,2, •■•,7), where 
T(q represents each component of the energy momen¬ 
tum tensor, {AA\ A, B, A^°\ ,G^ S \ F}. Since G^ 

vanishes identically outside the source, the homogeneous 
equation for the modified master variable should be un¬ 
altered by the transformation, 
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(e) c-» (e) c+E c ‘ G (w ( 5 °) 

2=1 


Now one may think that the factor \C\ 2 from the ex¬ 
pression for can be eliminated by using this degree 
of freedom of modifying the master variable. However, 
we will prove below that it is impossible. 

As a result of the transformation above, is mod¬ 
ified as 


A;(T) x(T) _ fc{M) 


Y. Ci ^co 
. 2=1 


(51) 


Since A'( M ) [( e )£] contains we cannot eliminate 

T(q jT . from unless d = 0 for i > 


4. For {A^ l \A,B}, one can use the conservation 
law cnu to eliminate T(i), r - Thus, the condition that 
\C\ 2 dT^\ ft 0 requires d ft 0 for i > 4, where 
« means the equality modulo |C| 2 . Then, we find that 
and G (s) arises in the expression for the modified 
only from J 2 i <3 a Tf t ^ r . Hence, the conditions 
for the coefficients of B a nd G^ to vanish modulo 
\C\ 2 become c\ ft -(wr 2 /\/2(l + A)(r — 2M))c3, and 
c 2 ft \/(1 + A)c 3 , respectively. Thus a possible modi¬ 
fication which might eliminate the factor |C|~ 2 from the 
expression for F T ) is restricted to 

AT (t) -a K^ - A (M) [/(r)(-(wr 2 /V2(l + A)(r - 2 M)) 

xA« + ^(T+A)A + H)], (52) 

with an arbitrary function /(?'). Then, a straightforward 
calculation shows that the factor |C|~ 2 cannot be elimi¬ 
nated by this transformation. Thus, the idea of introduc¬ 
ing a new master variable for even parity modes satisfying 
the RW equation does not work well for the purpose of 
metric reconstruction in time domain. 


IV. DISCUSSION 

In this paper we have introduced new master variables 
for the odd and the even parity cases. We call them, re¬ 
spectively, the modified Regge-Wheeler and Zerilli vari¬ 
ables. These variables satisfy the same Regge-Wheeler 
or Zerilli equation except for the source terms, which are 
composed of the matter energy momentum tensor. We 
have given the explicit expressions for the source term. 
The metric perturbation in the RW gauge is expressed 




































in terms of the master variables and the matter energy 
momentum tensor. The explicit formulas for the metric 
reconstruction were also written down. The important 
aspect of these modified variables lies in the fact that 
the frequency ui does not appear in the denominator in all 
the formulas to obtain the metric perturbation. Hence, 
there is no time integration except for the step solving 
the master equation. The most crucial point will be that 
/t' T ), the contribution to the reconstructed metric pertur¬ 
bation from the matter energy-momentum tensor, does 
not have w in the denominator. Therefore, the perturbed 
metric around a field point (t, r) is solely written in terms 
of the master variables if the energy-momentum tensor 
vanishes in the vicinity of the spherical shell containing 
this field point. This fact will be useful in the program 
to calculate the regularized self-force acting on a particle 
orbiting in the Schwarzschild spacetime. 

As mentioned earlier, in the Introduction, the full met¬ 
ric perturbation contains a divergent piece near the par¬ 
ticle location. To obtain a sensible expression for the 
self-force, we need to subtract the so-called “direct” part 
from the full metric perturbation before evaluating the 
expression of the force. But here is the “gauge prob¬ 
lem” . The full metric perturbation is obtained in Regge- 
Wheeler gauge but the “direct” part is evaluated in the 
harmonic gauge associated with the particle trajectory. 

Here we would like to propose an insight towards an 
alternative method to handle this gauge issue in the case 
of the Schwarzschild background. The basic idea is in¬ 
spired by the notion brought by Barack and Ori 0 . They 
stressed that the trajectory in the perturbed spacetime is 
gauge invariant although the expression for the self-force 
depends on the choice of gauge. On the other hand, the 
metric perturbations reconstructed from this gauge in¬ 
variant master variables depend on the choice of gauge, 
but the concepts of the perturbed geometry and hence of 
the geodesic on it are gauge invariant. Hence, naturally 
one may expect that the subtraction at the level of gauge 
invariant master variables is possible. 

A sketch of the new method is the following. The “di¬ 
rect” part of the metric perturbation h < - s ' ) can be calcu¬ 
lated in the harmonic gauge. We can use the recent ob¬ 
servation by Detweiler and Whiting E3 that can be 
modified so that it satisfies the Einstein equations. Since 
the method for the harmonic decomposition of the di¬ 
rect part is established by Mino, Nakano and Sasaki Q, 
the projection of this direct part to the gauge invariant 
master variable is possible by using the formulas 
m and iPlt . On the other hand, solving the RW equa¬ 
tion, we can directly calculate the master variable corre¬ 
sponding to the full metric perturbation, £(f ull '>. Then 
we subtract the direct part from to obtain the 

master variable that corresponds to the tail part, which 
we denote by ^ tal1 ). Since both and satisfy 

the RW equation with the same source, their difference 
Q tall '> satisfies the homogeneous RW equation. Hence, 
we can reconstruct the metric perturbation correspond¬ 
ing to the tail part from this regularized master variable 


£(*“**) by applying the formulas h^ M \ At this step the 
choice of gauge is unimportant as is explained in the pa¬ 
per by Barack and Ori 0 . Since the subtraction of the 
divergent part is done at the level of the gauge invari¬ 
ant variables, we would like to call this scheme the gauge 
invariant regularization. 

In the new scheme, using the variables introduced in 
this paper, the part depending on the master variable 
in the metric reconstruction formulas, hS M \ does not 
have uj in denominator as well as h ^ . Hence, when we 
know the behavior of the master variable corresponding 
to a homogeneous solution of metric perturbations in the 
vicinity of a spherical shell, we can reproduce the met¬ 
ric perturbations from the master variable there. If 
contained u> in denominator, the local information of the 
master variables near the shell were not sufficient to re¬ 
produce the metric perturbations. Therefore, the use of 
the new variables introduced in this paper is crucial for 
the gauge invariant regularization. 

This scheme still has a subtle point which requires fur¬ 
ther investigation. The method for the reconstruction 
of the metric perturbation does apply only for a solu¬ 
tion of vacuum Einstein equations. However, in the ac¬ 
tual computation, the direct part h^ is calculated in a 
power series expansion with respect to the separation £ 
between the source point and the field point, and this ex¬ 
pansion must be truncated at a certain order of £. Then, 
the truncated direct part does not satisfy the Einstein 
equations in general. Hence, we need a new invention 
to bypass this difficulty in order to realize this attractive 
idea of the gauge invariant regularization. We would like 
to return to this challenging issue in a future publication. 

In Sec. Ill we discussed the possibility of using a mas¬ 
ter variable for even parity modes which has the same 
potential for the master equation as in the case of odd 
parity modes. Such a variable is obtained by using the 
Chandrasekhar transformation. We wrote down the ex¬ 
plicit definition of this master variable in terms of the 
metric components, the master equation with the source 
terms and the metric reconstruction formulas. We found 
that the metric reconstruction formulas necessarily con¬ 
tain the Starobinsky constant including u> in the denom¬ 
inator. Therefore, the use of the even parity master vari¬ 
able that has the same homogeneous master equation as 
in the odd parity case unfortunately turned out not to 
be advantageous. However, the master variables consid¬ 
ered here are limited to those which are related via Chan¬ 
drasekhar transformation. We expect an even wider class 
of transformations in which we might find a more suitable 
variable for the purpose of metric reconstruction. 
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